Skip to main content

Full text of "A Detection of Gas Associated with the M 31 Stellar Stream"

See other formats

arXiv: 1506.00635vl [astro-ph.GA] 1 Jun 2015 

Accepted for publication in the Astrophysical Journal; May 28, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


Andreas Koch^, Charles W. Danforth^, R. Michael Rich®, Rodrigo Ibata'^, Brian A. Keeney® 
Accepted for publication in the Astrophysical Journal; May 28, 2015 


Detailed studies of stellar populations in the halos of the Milky Way and the Andromeda (M31) 
galaxies have shown increasing numbers of tidal streams and dwarf galaxies, attesting to a complicated 
and on-going process of hierarchical structure formation. The most prominent feature in the halo of 
M31 is the Giant Stellar Stream, a structure in extent along the sky, which is close to, but 

not coincident with the galaxy’s minor axis. The stars that make up this stream are kinematically 
and chemically distinct from the other stars in the halo. Here, we present HST/COS high-resolution 
ultraviolet absorption spectra of three Active Galactic Nuclei sight lines which probe the M3I halo, 
including one that samples gas in the main southwestern portion of the Giant Stream. We see two 
clear absorption components in many metal species at velocities typical of the MSI halo and a third, 
blue-shifted component which arises in the stream. Photoionization modeling of the column density 
ratios in the different components shows gas in an ionization state typical of that seen in other galaxy 
halo environments and suggests solar to slightly super-solar metallicity, consistent with previous 
findings from stellar spectroscopy. 

Subject headings: galaxies: halos — quasars: absorption lines — galaxies: individual (MSI) — galax¬ 
ies: ISM — galaxies: kinematics and dynamics 


The assembly of galaxies via hierarchical accretion - 
wherein smaller structures merge to create larger struc¬ 
tures - is now a well-established doctrine of galaxy 
formation, which only requires the action of gravity. 
This universal prediction of co smological simulations 
(e.g., iBovlan-Kolchin et al.ll20I^ is also well observed 
in the local universe, with simulations predicting a mul- 
titude of substructure s in the Milky Way (MW) halo 
(|Helmi fc White! [19991 1. Streams of stars and gas have 
been observed in the MW and have b een linked to 
previ o us accretion of s a tellite galaxies (iSearle fc Zinnl 
119781 : llbata et ^ 11994 iBelokurov et al.l 1200711 . The 
Magellanic Stream ( Mathewson et al.l 11974) observed 
in 21-cm emission is the result of the interaction be¬ 
tween the Large and Small Magellanic Clouds (LMC 
and SMC, respectively) and the MW and may, gi- 
gayears hence, result in the merger of the Magellanic 
and MW galaxies. On the s tellar side, the a dvent of 
the Sloan Digital Sky Survey (|York et al.l 1200011 has led 
to an explosion in the discovery of s t ellar streams in 
the MW halo (e.g.. lYannv et al.l[20031 IBelokurov et ahl 

2007HGrillmair 1 120091: iNewberg. Yannv. fc WillettI 120(1 


Martin et al.ll2013h . all of which presumably have their 

origin in minor mergers in the Galaxy’s past. 

t Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated by the Associa¬ 
tion of Universities for Research in Astronomy, Inc. under NASA 
contract NAS 5-26555. 

^ Zentrum fiir Astronomic der Universitat Heidelberg, Lan- 
dessternwarte, Konigstuhl 12, 69117 Heidelberg, Germany 
^ CASA, Department of Astrophysical and Planetary Sciences, 
University of Colorado, 389-UCB, Boulder, CO 80309 

® Physics and Astronomy Building, 430 Portola Plaza, Box 
951547, Department of Physics and Astronomy University of 
California, Los Angeles, CA 90095 

^ Observatoire de Strasbourg, 11, rue de I’Universite, F-67000, 
Strasbourg, France 

In the past decade our knowledge of the structure of 
our nearest large galactic neighbor, M3I, has under¬ 
gone a similar tran s forma tion. Deep, wide-held imag¬ 
ing by llbata et ^ (|2001ll revealed a large stream of 
stars extending from near the center of M3I several 
degrees to the south east close to the minor axis of 
the galaxy. Subsequent studies ((Eerguson et al.l 120021 : 
iMcConnachie et al.ll2003l l found that this Giant Stellar 
Stream (GSS) exists as a coherent structure more than 
4° on the sky, inclined at ^60° with respect to the line of 
sight, and extending over > 140 kpc. The GSS is aligned 
with the major M31 satellites M32 and NGC205 and, 
while it is tempting to associate the stream with those 
companions, modeling shows that it most likely repre¬ 
sents the trailing arm left from a recent, radial merger 
of a galaxy of 3 — 5 x 10® M 0 (Mori & Rich 2008, Eardal 
et al. 2013). Alternatively it may be the fossil remnant 
of a cannibalized satellite or a tidal bridge connecting to 
it s neighboring spiral M 33. 

iBrown et al.l ()2006l l found from Hubble Space Tele¬ 
scope (HST) color-magnitude diagrams that the Stream 
and halo helds are indistinguishable, which suggests that 
only one progenitor was responsible for the debris de¬ 
posited in inner halo regions. However, if dynamical 
mixing were efficient, it could erase the signatures from 
different infall sources. Further deep imaging and stel¬ 
lar spectroscopic studies have shown that the halo of the 
M 31-M 33 system is far more complicated than originally 
thought and hint s at an eventful formation history (e.g., 
llbata et ^1120071 : IMcConnachie et al.ll2009l l. In addition 
to the GSS, many other structures have been discovered 
within the halo of M31. Some of these structures are 
associated with the large a nd growing population o f An¬ 
dromeda satellite galaxies (iRichar dson et ^ 1201111 and 
globular clusters (|Mackev et al.l 1201011 . while others are 
less clearly connected. It has been suggested that many 
of the stream-like features in the halo could result from 


A. Koch et al. 

the heating of M31’s thin disk that has been ejected dur¬ 
ing a recent merger (e.g., Richardson et al. 2008; Bernard 
et al. 2015). 

We know relatively little about the progenitor of the 
GSS. Modeling of the encounter leads to mass estimates 
ranging from 3 x 10® Mq (Hammer et al. 2010) to 
< 5 X 10® Mq (Mori & Rich 2008); the most recent esti¬ 
mate of 3 X 10® Mq comes from Fardal et al. (2013), plac¬ 
ing the system at the mass of the LMC. It has been sug¬ 
gested that the lack of a compact remnant favors a disk¬ 
like progenitor. Furthermore, Fardal et al. (2008) argue 
that asymmetries in the surface brightness distribution 
across the GSS and its morphology in terms of a more 
compact core with an extended envelope can best be ex¬ 
plained if this progenitor disk galaxy was rotating. Based 
on the metal-rich nature and the age-metallicity relation 
of several halo fields with stream-like color magnitude 
diagrams Bernard et al. (2015) suggest an early-type 
galaxy as a candidate, which experienced rapid chemical 
enrichment to Solar abundances already hy z = 1. 

However, no neutral gas analogous to the Magellanic 
Stream is evident in the immediate vicinity of the GSS 
(Braun et al. 2003). A recent re-analysis of the Hi 
data by Lewis et al. (2013) showed the presence of gas 
filaments, bridging M31 and the near-by spiral M33, 
although they are significantly offset from the stellar 
stream. Lewis et al. (2013) also identified a number of 
individual gas clouds, but noted a remarkable lack of any 
spatial correlation between stellar and gaseous substruc¬ 
tures, arguing for several dynamical processes at play 
that affect both components in different ways. The only 
gas associated with the GSS thus appears at the point 
where it passes the disk of M31. 

Absorption-line spectroscopy provides a different and 
complementary tool for analyzing the structure and kine¬ 
matics of circumgalactic matter (e.g., Rao et al. 2013). 
Stellar spectroscopy provides a representative view of 
stellar kinematics and overall metallicities, while absorp¬ 
tion spectroscopy of background sources gives kinematic 
and chemical information about the gas phase in the sys¬ 
tem. In this paper, we present far-ultraviolet (FUV) 
spectra of three active galactic nuclei (AGN) obtained 
with the Gosmic Origins Spectrograph (COS) onboard 
the HST. These three sight lines probe the M31 halo at 
projected distances of 26-51 kpc in three different quad¬ 
rants of the system (Fig. 1). 

We see strong absorption from gas near the systemic 
velocity of M 31 in all three sight lines. The sight line that 
probes the GSS additionally shows clear absorption at 
more negative velocities, near the systemic velocity of the 
stellar members of the Giant Stream. While Lehner et 
al. (2015) associate such absorbers with the hot, gaseous 
corona of M31, no connection to the stellar component 
of the GSS has been made. 

Our measurements of the M31 and GSS-velocity gas 
are presented in §3. In §4, we compare the observed gas- 
phase detections with the stellar kinematics in the same 
field and a kinematic model from the literature. Here, 
we also analyze the absorbers by application of basic pho¬ 
toionization models. Our conclusions are presented in §5. 


We base our analysi s on HST/COS (jGreen et al.ll2012l : 
lOsterman et al.l|2011l) observations of three AGN sight 


50 0 -50 


-0.60 < [Fe/H] < 0.00 

f (degrees) 

Figure 1. The locations of the three AGN used in this study 
overlaid on a density map of metal-rich stars (after Ibata et al. 
2014). Notice that Q0045-|-3926 lies near the edge of the Giant 
Stellar Stream. Here, we have masked out the inner disk of M31 
within a radius of 17 kpc. The inner halo out to a radius of 30 kpc 
and at an axis ratio of 3:5 is displayed at lower contrast to render 
visible the inner halo substructure. 

lines, which pass near M31 (Fig. 1). Of the three sight 
lines, the one near the GSS (Q 0045-1-3926) is notewor¬ 
thy as it lies 0.2° (~2.7 kpc) East of the main stellar 
body of the GSS. The predominance of stellar stream 
features and complexity in all regions of the halo is well 
documented (e.g. McGonnachie et al. 2009). Three 
Seyfert galaxies were observed in 2009, 6-7 October as 
part of Program 11632 (P.I. R.M. Rich). Each observa¬ 
tion consisted of a single exposure in the GOS/G130M 
(1135 < A < 1450 A; 6.00 ksec) grating and one to three 
exposures with GOS/G160M (1400 < A < 1795 A). In 
addition, Q 0045-1-3926 was observed in late 2010 (P.I. N. 
Arav). Observation details are given in Table 1. 

The exposures were obtained from the Mikulski 
Archive for Space Telescopes. The calibrated, one¬ 
dimensional spectra for each target were next coadded 
with the custom IDL procedures described in detail by 
iDanforth et al.l (|2010f) . The coaddition routines auto¬ 
matically scales exposures taken during different epochs 
to take into account source variability and stitch G130M 
and G160M data together into a continuous spectrum 
over the range of 1150-1750 A. Since the two objects 
HS 0033-1-4300 and HS 0058-1-4213 were observed with a 
single GOS grating position, each dataset features a pair 
of ^ 20 A gaps in the spectral coverage, corresponding 
to the gap between the two far-UV detectors. 

We ht continua to each of the data sets using a semi- 
automated line identification and spline-fitting technique 
as follows. First, the spectra are split into 5-10 A seg¬ 
ments. Gontinuum pixels within each segment are iden¬ 
tified as those for which the signal-to-noise (S/N) ratio 


Gas in the M 31 Stellar Stream 


Table 1 

Summary of COS Observations 

Target Name 

Alternate Name 






Pm 31 
[kpc] ^ 

Date of 
of observation 





Q 0045+3926 

lO And 
















HS 0033+4300 










HS 0058+4213 

J 0101+4229 










^ Minimum impact parameter from sight line to center of M31 based on the scale 228 pc arcmin ^ and mean 
distance 792 di 193 kpc. 

^ Program 11632, P.l. Rich 
^ Program 11686, P.L Arav 

vector is less than 1.5 a below the median S/N value 
for all the pixels in the segment. Thus, absorption lines 
(flux lower than the segment average) are excluded, as 
are regions of increased noise (error higher than segment 
average). The process is iterated to convergence and a 
spline is fitted through the mean continuum pixel value in 
each segment. We check the continuum fit of each entire 
spectrum manually and the continuum region identifica¬ 
tions are adjusted as needed. The continuum identifica¬ 
tion and spline-fitting processes work reasonably well for 
smoothly varying data, but they were augmented with 
piecewise-continuous Legendre polynomial fits in a few 
cases. In particular, spline fits perform poorly in re¬ 
gions of sharp spectral curvature, such as the Galactic 
Lya absorption and at the peaks of cuspy emission lines. 
Complete details on the process are given in Danforth et 
al. (2014). 

The region around M 31 has previously been mapped 
in the 21 cm emiss ion line of H I by several authors. 
iThilker et al.l (|2004ll conducted a Gree n Bank Telescope 
GBT ) survey of the region, while iBraun &: ThilkeJ 
2004fl present a deep survey in the M31 region made 
with the Westerbork Synthesis Ra dio Telescope (WSR T) 
operated in total-power mode (see IBraun et al.ll2003l for 
details). While the absorption-line data are much more 
sensitive than emission, the COS data give no informa¬ 
tion on the neutral hydrogen column density since the 
MW damped Lya absorption profile extends over thou¬ 
sands of km s“^. 

Unfortunately, the velocity calibration of the UV spec¬ 
tra is uncertain to a much larger degree than radio emis¬ 
sion data. HST spectra are recorded in a heliocentric 
velocity frame with an on-board wavelength calibration 
lamp. However, the precise wavelength solution varies by 
^0.2 A (^2 resolution elements) over the detector and 
the zero point is uncertain by ^^15 kms“^ due to inac¬ 
curacies in centering the target in the COS aperture in¬ 
herent to the target acquisition procedure. The coalition 
process cross-correlates strong absorption lines in each 
detector segment in each exposure with each other be¬ 
fore coaddition, so the relative wavelength solution over 
a few Angstroms is reasonably accurate. 

However, no cross-segment wavelength correlation is 
attempted. We correct for these effects to first order 
by aligning narrow low-ionization metal absorption lines 
in the data in velocity space to the H I 21-cm profiles 

from the Leiden-Argentine-Bonn survey (|Kalberla et al.l 
[200^ . The H I emission toward all three sight lines 
shows a peak at 0 < ulsr < +3 km s“^ with significant, 
distributed emission toward negative velocities. Fits to 
the Q 0045-1-3926 COS data show the reddest absorp¬ 
tion component in several ions at ^ —20 km s“^. Poor 
data quality in the other two spectra hampers this analy¬ 
sis rather severely, but the strongest Galactic absorption 
appears to be at —40 ^ v < —30 km s“^. These values do 
not vary considerably with rest wavelength, so we adopt 
a uniform velocity shift of -1-20 km s“^ for Q 0045-1-3926 
and -k35 km s”! for HS 0033-^4300 and HS 0058-^4213 
to bring the data approximately into the Local Standard 
of Rest (LSR) reference frame. Overall, we estimate a 
systematic velocity uncertainty of ±15 km s“^ inherent 
in the absolute velocity of any absorption measurements 

We measure absorption features in three ways. 
Firstly, unsaturated absorption features with clear com¬ 
ponent structure can be htted via the apparent optical 
dept h (AOD) and prohle-weighted mean quantity meth - 
ods (|Savage fc SeinbadilllOOltTSembach fc Savaggll992ll . 

Secondly, multi-component Voigt profile fits often pro¬ 
vide more accurate measurements in cases of mildly sat¬ 
urated lines. The significance of any individual line fit 
is calculated based on the local S/N ratio, the Doppler 
b-parame ter, and fitted equiv alent width via the algo¬ 
rithms of iKeenev et al.l (|2012ll . Finally, we make use of 
simultaneous multi-component Voigt profile fits to mul¬ 
tiple transitions of the same species, e.g., Si II 1260.42, 
1193.29, 1190.42, and 1526.70 A. These techniques pro¬ 
duce equivalent results for absorption profiles of mod¬ 
erate strength, but each has its advantages in different 
circumstances: Voigt fitting, which is generally preferred, 
can accurately measure lines which are moderately sat¬ 
urated (r ~2-3). The down-side of this technique is 
that a velocity component structure must be assumed. 
AOD measurements are well-suited for shallow, noisy ab¬ 
sorbers or where the component structure is uncertain, 
while this method tends to under-estimate column den¬ 
sities for very strong absorbers. The latter does not pose 
a concern in our only mildly saturated absorbers. For 
details on the consistencies between the methods, we re¬ 
fer the reader to Danforth & Shull (2008; and references 
therein). Measurements from all three techniques are 
performed on the data and the resulting values are pre- 


A. Koch et al. 

sented in Tables 2 and 3, where we label the absorber 
“components “ by the error-weighted mean of all species. 

3.1. Q 0045+3926 

The spectrum of Q 0045-1-3926 shows five distinct ab¬ 
sorption components in the range —400 < vur ^ 
0 km s“^ in many metal absorption lines from low- 
ionization Al II and Fe II to highly-ionized C IV tran¬ 
sitions (Fig. [5]). The strong absorption at n « 0 kms“^ 

-600 -400 -200 0 200 

Velocity (LSR, km/s) 

Figure 2. Normalized absorption profiles in various transitions 
along the Q 0045+3926 sight line. Data are binned to three pixels 
(~40% of the COS resolution) in all cases. The (small) normalized 
error vector is shown in gray. Vertical dashes correspond approxi¬ 
mately to the velocity component centroids discussed in the text. 

is similar to that seen in all other extragalactic UV ab¬ 
sorption lines and can be unambiguously associated with 
absorption from the MW disk gas. Other absorption, 
often blended with the MW gas, appears at —70 
km s“^. This component could plausibly be associated 
with the M31 halo given the overlap of MW and M31 
stellar velocities discussed in Sect. 4.1 below. However, 

this absorption is also consistent with infalling or out¬ 
flowing gas within the Milky Way halo in the form of 
the high velocity clouds (HVCs) and we deem this to be 
a more likely explanation (e.g., Wakker & van Woerden 
1997; Lehner et al. 2012). 

At n = —170 and —230 km s“^, two strong absorp¬ 
tion components can be seen in transitions of Si II (1193, 
1190, 1526 A), Si III (1206.5 A), Si IV (1393, 1402 A), 
and C IV (1548, 1550 A). The strong Si II 1260 A line 
is blended with weaker S II 1259 A absorption from the 
MW at that velocity, though the observed profile is con¬ 
sistent with the same absorption seen in three weaker 
transitions of the same ion. The strong C II 1334.5 A 
line is similarly complicated by HVC absorption from 
the excited-state C II* 1335.7 A transition, but is again 
consistent with the absorption component structure seen 
at u ~ —230 km s“^ in other metal ions. Notably, the 
relative strengths of the two lines in low ions as com¬ 
pared with higher ions in the same species (e.g., Si II 
and Si IV) show opposite trends implying that the ma¬ 
terial at n ~ —170 km s“^ is less ionized than the gas at 
V —230 km s“^. 

Prominent absorption is seen at n ~ —370 km s“^ in 

Si II 1260 A, C II, and Si III and weakly in other Si II 
transitions as well as both Si IV lines and C IV 1548 A. 
This is a strong, distinct, system that is ^ —70 kms“^ 
relative to the M31 systemic velocity, and we associate 
it with the GSS, despite the lack of obvious H I emission 
in the immediate GSS, given the lower column densities 
that our observations sample. We will further discuss its 
possible origin and connections to stellar kinematics in 
its surroundings in Sect. 4.1. 

We note that all the detected three high velocity com¬ 
ponents could in principle arise in Milky Way HVGs as 
well, although we deem this very unlikely: MW halo 
gas clouds are typically observed at velocities up to 
^ 200 km s“^, and the ex treme wing of the distribution 
reaches up to 400 km s“^ (|Wakker fc van WoerdenifIMI: 
iPutman et al.l 1^00211 . Some weak metal-line absorption 
systems have been detected in the MW halo that are not 
associated with 21-cm emission. These detections have 
come in both low-ioniza t ion sp ecies such as Ga II and 
Na I (iRichter et al.ll2005l l2009li and higher ions, such as 
G IV and O VI ( Sembach et al.l[^003ll . 

3.2. HS 0033+4300 and HS 0058+4213 

The data for HS0033-f4300 and HS0058-f4213 are 
considerably noisier than those of Q 0045-1-3926 due to 
fainter source fluxes and shorter exposure times. De¬ 
tailed hts to the absorption profiles in these two sight 
lines are difficult, but similar trends are qualitatively 
seen. Fig. [3] shows absorption aX v -1-15 and —70 
km s“^ in all bands which we ascribe to MW disk and 
halo gas. Both sight lines show absorption at n ^ —180 
to —170 km s“^ in the moderately and highly-ionized 
species as well, presumably tracing M 31-velocity gas de¬ 
spite their considerable distance from the main body of 
the disk. Both the MW and M31 absorption may plau¬ 
sibly be sub-divided into additional components as seen 
in the higher-quality data of Q 0045-1-3926. 

HS 0033-1-4300 shows a low-significance absorption fea¬ 
ture in several ions (Si III, G IV) consistent with v ^ 

Gas in the M 31 Stellar Stream 


Table 2 

Species Measurements Toward Q 0045+3926 





log N 


V ^ —170 km s ^ component 


-170 ±5 

20 : 

12.7 I*" 

< Ifj, AOD, Voigt 


-158 ±4 

21 ± 7 

13.39 ±0.06 

A1190, 1193; AOD, Voigt 


-162 ±2 

19 ±4 

13.17 ±0.11 

AOD, Voigt 

Si IV 

-160 ±7 

22 ±5 

12.92 ±0.12 

doublet; AOD, Voigt 


-130 : 

30 ; 

14 : 

Voigt only, strongly blended 


-173 ±2 

18 ±3 

13.57 ±0.03 

AOD (doublet), Voigt (A1548 only) 


< 12.97 

doublet; AOD upper limit 

V ^ —230 km s ^ component 


-236 ± 7 

20 : 

13.1 I*" 

AOD, Voigt la) 


-222 ± 16 

24 ±7 

12.99 ±0.29 

A1190, 1193; AOD, Voigt 


-228 ±4 

30 ±5 

13.16 ±0.04 

AOD, Voigt 

Si IV 

-225 ±4 

28 ±4 

13.22 ±0.10 

doublet; AOD, Voigt 


-227 ±3 

24 ± 5 

13.89 ±0.02 

AOD, Voigt 


-230 ± 3 

25 ±7 

13.98 ±0.07 

doublet; AOD, Voigt 


-230 ± 5 

10 : 

12.97 ±0.12 

~ 3(7, doublet; AOD, Voigt 

V —370 km s ^ component 


-374 ± 7 

20 : 

13.4 I*" 

AOD, Voigt 2a) 

Si II 

-359 ±4 

24 ±2 

12.88 ±0.10 

A1190, 1193,1260; AOD, Voigt 


-355 ± 2 

35 ±6 

13.09 ±0.04 

AOD, Voigt 

Si IV 

-390 ± 9 

40 : 

12.62 ±0.11 

A1393 only; AOD, Voigt 


-351 ± 5 

29 ±6 

14.12 ±0.05 

AOD, Voigt 


-378 ± 1 

34 ±5 

13.20 ±0.01 

AOD (doublet), Voigt (A1548 only) 


-393 ± 16 

30 : 

12.65 ±0.06 

< 3a; A1238 only; AOD, Voigt 

^ All line profiles have been shifted by +20 km s ^ as discussed in the text. 

^ The 3a detection limit for 01 lines is log N < 13.54. Marginal, lower-significance mea¬ 
surements are presented for completeness. 

Table 3 

Species Measurements Toward HS 0033+4300 and HS 0058+4213 





log N 


[cm J 

HS0033±4300, v 

^ —170 km s 

^ component 


< 14.15 

upper limit 


< 13.2 

A1190,1193,1526; simulfit 


-142 ± 3 

68 ±4 

13.5 ±0.3 

AOD, Voigt 

Si IV 

-177 ±2 

53 ±18 

13.5 ± 0.1 

A1402 only; AOD, Voigt 


-155 : 

65 : 

14.5 : 

very uncertain, Voigt, AOD 


-176 ± 4 

35 ± 7 

14.13 ±0.07 

doublet; AOD, Voigt 


< 13.56 

upper limit 

HS0033±4300, v 

~ —350 km s“ 

^ component 


< 14.15 

upper limit 


< 13.2 

upper limit 


-350 ± 5 

22 : 

12.9 : 

< 3(t Voigt, AOD 

Si IV 

A1402; instrumental feature 


-359 ± 11 

> 50 

14.3 ±0.1 

AOD, Voigt 


-323 ±21 

75 ±12 

13.77 ±0.09 

doublet; AOD, Voigt 


< 13.56 

upper limit 

HS0058±4213, v 

^ —180 km s“ 

^ component 


< 13.9 

upper limit 


-188 ± 11 

37± 15 

13.8 ±0.2 

A1190, 1193,1526; AOD, Voigt 


-179 ± 10 

15 ±5 

14.5 ± 0.3 


Si IV 

-183 ± 8 

38 ±7 

13.45 ±0.18 

doublet; AOD, Voigt 


-183 ±2 

33 ±6 

14.5 ± 0.1 

AOD, Voigt 


-162 ± 8 

53 ±14 

14.15 ± 0.13 

doublet; AOD, Voigt, simulfit 


< 13.4 

upper limit 

^ All line profiles have been shifted by +35 km s ^ as discussed in the text. 


A. Koch et al. 

-600 -400 -200 0 200 -600 -400 -200 0 200 

Velocity (LSR, km/s) Velocity (LSR, km/s) 

Figure 3. Same as Fig. [2] for HS 0033+4300 (left) and HS 0058+4213 (right). Data are binned by three pixels. The stronger Si IV line 
(A = 1393) is blended in the left panel with strong intrinsic absorption, so the weaker A = 1402 line is shown. Vertical lines are the five 
components seen in Q 0045+3926. 

Gas in the M 31 Stellar Stream 


—350 km gas, though there is little evidence for 
this absorption in lower-ionization species. This con¬ 
trasts Q 0045-1-3926 where the strongest absorption at 
this velocity is seen in the low-ionization species Si II, 
Si III, and C II. No convincing absorption is seen in the 
HS 0058-1-4213 sight line at n < —300 km s“^. 


Our examination of the three absorption spectra finds 
that all have in common Galactic components, and com¬ 
ponents at around —170 to —230 kms“^ that are plau¬ 
sibly associated with the M31 disk. 

4.1. Gas Kinematics and Comparison with Stellar 

Q 0045-1-3926 stands out in exhibiting a unique com¬ 
ponent at —370 kms“^ not seen in the other two spectra, 
although there is a low-significance absorption at —350 
kms“^ toward HS 0033-1-4300. The simulations of Mori 

6 Rich (2008) suggest that a stream component would 
be found in this NW region, but the putative components 
appear too weak for a firm association. As the sightline 
towards Q 0045-1-3926 passes ^2.7 kpc from the GSS 
we now turn to the question of whether observed stel¬ 
lar radial velocities or modeled stellar stream velocities 
are consistent with the —370 kms“^ component in the 
gas. Fig. 4 (bottom panel) shows stellar radial veloci¬ 
ties derived from Keck/DEIMOS measurements within 

7 kpc of Q0045-1-3926 (Ibata et al. 2015 in prep.); the 
stream shows a clear signature with a peak near —500 
kms“^, well separated from the gas component. Such a 

Velocity v [km s 

Figure 4. Observed stellar radial velocity distributions in the 
immediate vicinity (7 kpc) of two of our AGN sightlines, where 
available: HS 0058-1-4213 (top panel) and Q 0045-1-3926 (bottom 
panel), taken from Ibata et al. (2004) and Ibata et al. (2015, in 
prep.). The bottom field contains a strong signature of the Giant 
Stream at —500 kms“^. Solid and dashed lines indicate the detec¬ 
tion of gas absorption at those velocities, as listed in Tables 2,3. 

difference is in line with the lack of any obvious spatial 

correlation between stellar and gas velocities reported 
by Lewis et al. (2013). The presence of a strong stellar 
velocity gradient within the GSS across its 4.5° extent 
has already been noted by Ibata et al. (2004), ranging 
from —300 kms“^ near its tip to —500 kms“^ towards 
the disk of M31 (see also Koch et al. 2008; Gilbert et al. 
2009). While the gas filaments also display a strong kine¬ 
matic gradient, there are two marked differences (Lewis 
et al. 2013): firstly, gas is only detected with a signif¬ 
icant parallel offset to the stellar stream and secondly, 
gas near the disk shows the systemic velocity of M31, at 
—300 kms“^, transiting to more positive velocities as it 
approaches M 33. Koch et al. (2008) and Gilbert et al. 
(2009) also identify a kinematic signature of the GSS at 
—440 kms“^ and thus much closer to the gas velocity 
of our present work. The respective fields (“a3”) are, 
however, located at ^30 kpc from Q 0045-1-3926, adding 
further complexity to the interpretation of velocity gra¬ 
dients and offsets between stream and stellar dynamics 
and spatial coincidences. 

The top panel of Fig. 5, in turn, shows the Mori & Rich 
(2008) model distribution of velocities for the M 31 halo 
and stream particles within 5 kpc around Q 0045-1-3926. 
The models had been tailored to reproduce the morphol¬ 
ogy and kinematics of the GSS and distinguished parti¬ 
cles from the disk/bulge, halo, and the accreted satellite. 
First, we note that no disk or bulge particles, ejected 
during the merger, are found in this region. Halo stars 
in the simulation show, per construction, the broad dis¬ 
tribution of stellar velocities that also include the three 
absorbers we observe in the gas phase and in the stel¬ 
lar samples. Turning to the satellite particles that make 
up the stream in the simulations, there is a clear coin¬ 
cidence of the broad distribution with the observed gas 
velocities at —370 km s“^. Moreover, the stream ve¬ 
locity distribution shows kurtosit0 toward less negative 
velocities with a peak at —300 kms“^. 

We now turn to the case of HS 0058-1-4213. Fig. 4 (top) 
shows that the respective stellar sample was dominated 
by MW foreground dwarfs above ^ —100 kms“^, which 
level off into the M31 extended disk and halo profile be¬ 
low ^ —150 kms“^; a tail of stars is found at —400 
kms“^ < V < —300 kms“^, which represents the canon¬ 
ical halo population. The observed gas absorption in this 
sight line at —180 kms“^ is also consistent with veloci¬ 
ties of individual gas clouds in the M31 halo (e.g., Lewis 
et al. 2013) so that both stellar and gas velocities in 
this region probe M31’s halo with no obvious connection 
to any substructures. We note that no Stream particles 
are found in this region of the Mori & Rich (2008) sim¬ 
ulations and the only dominant component is the broad, 
halo model peak at —300 kms“^. 

Finally, while no stellar velocities are available in the 
vicinity of HS 0033-1-4300, the simulations predict that a 
rather diffuse stream component should be seen at highly 
negative velocities and another peak that blends into the 
MW foreground, above ^0 kms“^ (Fig- 5, bottom). The 
gas absorption seen at —170 and, more weakly, —350 
kms“^ thus appears to be consistent with typical halo 

® Prom the 140 stream particles within 5 kpc of Q 0045-1-3926, 
we find a mean velocity of —340 kms“^ with a Icr-dispersion of 40 
kms“^, skewness of 1.2, and kurtosis of 3.1. 

A. Koch et al. 

o 0.1 


I 0.05 


0) U 


1 0.06 





-800 -600 -400 -200 0 200 

Velocity v [km s“^] 

Figure 5. Same as Fig.[^ but fo r the r adial velocity distributions 
from the models of IMorifc Richl ()2QQ8l h The gray-shaded and 
clear histograms are drawn from the halo, and stream particles, 
respectively, while vertical lines indicate the detection of gas ab¬ 
sorption at those velocities in the respective sight lines, as listed 
in Table 2. No stream particles are present in the vicinity of HS 

All three sight lines thus show absorption at near 
the systemic velocity of the M31 halo (v ^ —170 
km s“^), considering the stellar halo’s broad velocity dis¬ 
persion (at ^85 kms“^; Reitzel & Guhathakurta 2002). 
Q 0045-1-3926 shows two clear components separated by 
^ 60 km s“^ and it is clear from the ratios of high 
ions to low ions that these components have different 
ionization properties in different sightlines, in that the 
V ~ —170 kms“^ system shows weaker CIV and stronger 
low-ionization lines than the one v ^ —230 kms“^. Sim¬ 
ilarly, towards HS 0033-1-4300, C IV and Si IV absorption 
is seen near v ^ —176 km s“^ while C II and Si III show 
absorption centroids at —150 km s“^. 

In conclusion, we note that stars are found with veloc¬ 
ities consistent with all the three absorbing systems in 
the vicinity of the AGN sight lines, which is a manifes¬ 
tation of the large velocity dispersion of the M31 halo 
system. Based on the observed stell ar velocities and the 
comparison with the simulations of IMori fc Richl (|2008ll 
we argue that the lowest-velocity component present in 
the AGN spectra are associated with the Stream itself. 

4.2. Modeling Physical Conditions 

With our measurements of many metal ion column 
densities, it is possible to constrain the physical condi¬ 
tions in the absorbing gas with simple photoionization 
models. To this end, we compare measured column den¬ 
sities with a grid of GLOUDY models dev eloped for cir- 
cumg alactic absorbers (described fully in iKeenev et al.l 
1201311 . Briefly, a slab of gas is illuminated by an ion¬ 
izing c ontinuum given b y the extragalactic radiation 
held of iHaardt fc Madaul (120121) . As in Keeney et al. 
(2013), we use a purely-extragalactic (AGN) background 
since all three of our sight lines probe regions outside 
M31’s proximity zone where ionizing hux from early- 





type stars is expected to have an effect (Uprox 20 kpc; 
iGiroux fc Shulll [1997I1 . We calculate a grid of models 
where free parameters ionization parameter U = n^/nn 
and metallicity Z are varied to provide a grid of model 
conditions. Best-Ht model parameters for the three ab¬ 
sorption components toward Q 0045-1-3926 are listed in 
Table 4. 

Measurements in three adjacent ionization states of a 
metal provide a strong constraint on the ionization pa¬ 
rameter. In these data the strongest constraints are pro¬ 
vided by Si II/III/IV, though GII/IV ratios are generally 

H I absorption at the velocity of M31 is unrecover- 
ably blended with Galactic absorption which makes the 
metallicity of the absorbing gas difficult to constrain di¬ 
rectly. We set an upper limit of Nui < 5 x 10^^ cm“^ 
since the Q 0045-1-3926 sight line lies out side the low¬ 
est H I 21 cm column density contour of iThilker et al.l 
(|2004f) . However, we can infer Vhi for each compo¬ 
nent by hnding the range of total column density Nh 
over which all of the silicon ion column densities in a 
given component can be explained simultaneously in a 
single, purely-photoionized phase of a given metallicity 
Z. This gives values for Nh and Z. These quantities, 
along with the hydrogen neutral fraction {fni, uniquely 
given by the ionization parameter above) determines the 
neutral hydrogen c olumn density Ah t all within fairly 
large uncertainties (IStocke et al.lfeOOTl) . The equilibrium 
temperature, T, is a value returned by the best-fit U, Z 
model. A line-of-sight cloud size is then inferred in turn 
as Del = riH Jhi /Nbi ■ Assuming a spherical cloud of 
constant density, we calculate a cloud mass Md- 

We note that many assumptions have gone into these 
models and quantities farther down the chain of logic 
should be taken correspondingly with appropriate cau¬ 
tion. Modeled ionization parameters are probably rea¬ 
sonably accurate while inferred cloud mass may be inac¬ 
curate by orders of magnitude or more. 

The GLOUDY model best fit parameters are log U k, 
—3.0 for both the GSS gas and the component at —170 
km s“^, typical of conditions seen in the MW halo (e.g., 
Shull et al. 2011). The gas at —230 km s“^ appears to be 
more highly ionized (log U ~ —2.6) and thus typical of 
gas in “highly-ionized” HVGs, which is in agreement with 
the qualitative estimates of the respective line strengths 
seen in Fig. 2. 

The derived metallicities for the three absorption com¬ 
ponents along the Q 0045-1-3926 sight line are all consis¬ 
tent with solar or slightly super-solar metallicity, though 
with large uncertainties. In particular, for the two M 31- 
velocity components which are kinematically unrelated 
to the GSS and presumably trace gas in the halo of M 31 
(see also Lehner et al. 2015), these values are higher 
than the average stella r metalliciti e s in t he, generally 
more metal-poor, halo. iKoch et ^ (1200811 found stellar 
metallicities of [Fe/H]ss —1.4 dex in fields at 20-40 kpc 
from the galactic center. Stars in the Giant Stream, how¬ 
ever, are known to be more metal-rich than the surround¬ 
ing halo stars (Brown et al. 2006; Koch et al. 2008). 
For instance, Ibata et al. (2001) reported on an average 
[Fe/H] of approximately —0.7 dex, while Bernard et al. 
(2015) list a median metallicity of —0.35 dex based on 
their HST photometry. Moreover, the metallicity distri- 

Gas in the M 31 Stellar Stream 


Table 4 

Derived Physical Quantities from CLOUDY Models 


log AhG 

log U 

log Z/Zq^ 

log riH 
[cm“^ ] 



log /hi 


[kpc], line of sight 



Q 0045-1-3926, -370 km s”! 

15.96 ±0.33 

QC + O.IS 







Q 0045-1-3926, -230 km s”! 

15.36 ±0.30 

o f?o+0.26 







Q 0045-1-3926, -170 km s”! 

16.02 ± 0.48 

— 2 








HS 0033-1-4300, -350 km s"! 

17.3 ±0.4 


_1_0 1+0-3 






HS 0033-1-4300, -170 km s”! 

15.9 ±0.9 

_1_0 2^^'^ 






HS 0058-1-4213, -180 km s”! 

17.4 ±0.3 

-3.1 ±0.1 

—0 7+0-3 
‘ -0.2 






A^hi is determined by assuming that Sill/III/IV all arise in a single photoionized phase. 

For this, CLOUDY adopts the Solar metallicity of log Zq =0.012 from Asplund et al. (2005). 

Mci is the cloud mass (in solar masses) assuming a spherical cloud with density riff and diameter D^i- 

butions of Gilbert et al. (2009) peak at —0.2 dex and in¬ 
dicate a significant population of stars with roughly Solar 
metallicity in fields dominated by GSS material. These 
findings support our measurement of Solar metallicity in 
the stream gas, especially if one considers all the dif¬ 
ferent, employed measurement techniques in the studies, 
from stellar calcium triplet spectroscopy, to color magni¬ 
tude diagram fitting, and our modeling of the gas phase 

Given the poor quality of the data for the 
HS 0058-1-4213 and HS 0033-1-4300 sight lines we hesitate 
to place much weight on their model results. However, 
they are presented in Table 4 for completeness. 


We have investigated three sight lines in the M31 halo 
using HST/GOS and AGN as background sources. All 
of the sight lines exhibit components that we associate 
with the Milky Way, and as well, a component at —170 
to —230 kms“^ that we assign to the disk of M31. In 
the sight line of Q 0045-1-3926, which is projected ^2.7 
kpc from M31’s Giant Stellar Stream, we find S II and 
C II absorption at —370 kms“^, the most negative de¬ 
tected velocities in our sample. Here, we argue that 
the —370 kms“^ component in Q 0045-1-3946 is certainly 
associated with structures in M31, most likely with the 

The bulk of stellar velocities in the GSS covers a broad 
range, reflecting a strong radial velocity gradient across 
its body. Stars in the immediate vicinity of the AGN 
sightline move at near —500 kms“^, while other stream 
components have been found at velocities closer to our 
observed gas velocity (Ibata et al. 2004, 2015; Koch et al. 
2008; Gilbert et al. 2009), albeit at larger projected dis¬ 
tances. Aided by simulations of the merger that formed 
the GSS, Fardal et al. (2008) suggested that many of 
its morphological features could best be reproduced with 
a rotating disk-progenitor. In that case, a detachment 
of the stellar from the gas velocities would occur nat¬ 
urally. It is obvious that different mechanisms were at 
play, contrasting the tidal disruption of the stellar com¬ 
ponent of the stream progenitor with the gas removal 
by ram-pressure stripping in M31’s hot circumgalactic 
corona (Mayer et al. 2006, Lewis et al. 2013; Lehner 
et al. 2015). Likewise, the remarkable lack of coherence 
between stellar structures and presently detected large- 
scale gas in the halo emphasizes the complexity of the 
GSS formation and dynamics (Lewis et al. 2013). 

Further evidence for a stream-gas connection comes 
from the putatively high metallicity we derived for the 
gas phase. The halo of M 31 has long been known to have 
stars with metallicity similar to that of 47 Tuc (Mould 
& Kristian 1983; Rich et al. 1996). Even higher metal- 
licities were found in larger surveys, allowing for a bet¬ 
ter characterization of the main halo polluters. In fact, 
the HST-based age-metallicity relation of Bernard et al. 
(2015) showed that halo fields rich in substructures ex¬ 
perienced very rapid chemical enrichment so that, by a 
redshift of z = I, Solar metallicities were in place. This 
kind of evolutionary history is typically found in early- 
type (dwarf) galaxies (e.g., Layden & Sarajedini 2000; 
Bonifacio et al. 2004; Hendricks et al. 2014), allowing 
one to place further constraints on the progenitor’s iden¬ 
tikit. This emphasises the importance to increase the 
density of observational sampling of absorbers within 
the inner 30 kpc of the M 31 halo. This will be eased in 
the future by the increasing number of known AGN at 
suitable redshifts and magnitudes behind the M31-M33 
system, as listed, e.g., by Huo et al. (2013). 

AK gratefully acknowledges the Deutsche Forschungs- 
gemeinschaft for funding from Emmy-Noether grant Ko 
4161/1. The authors would like to thank K.R. Sembach 
for designing the observing program and D.B. Reitzel for 
helpful discussions. We also thank the anonymous referee 
for a swift and constructive report. 


Asplund, M., Grevesse, N., & Sauval, A. J. 2005, Cosmic 
Abundances as Records of Stellar Evolution and 
Nucleosynthesis, 336, 25 

Belokurov, V., Evans, N.W., Irwin, M.J., et al. 2007, ApJ, 658, 

Bernard, E. J., Ferguson, A. M. N., Richardson, J. C., et al. 2015, 
MNRAS, 446, 2789 

Bonifacio, P., Sbordone, L., Marconi, G., Pasquini, L., & Hill, V. 
2004, A&A, 414, 503 

Boylan-Kolchin, M., Springel, V., White, S. D. M. & Jenkins, A. 
2010, MNRAS, 406, 896 

Braun, R., Thilker, D., & Walterbos, R. A. M. 2003, A&A, 406, 

Braun, R., & Thilker, D. A. 2004, A&A, 417, 421 
Brown, T. M., Smith, E., Ferguson, H. C., et al. 2006, ApJ, 652, 

Danforth, C. W., & Shull, J. M. 2008, ApJ, 679, 194 
Danforth, C. W., Keeney, B. A., Stocke, J. T., Shull, J. M., & 
Yao, Y. 2010, ApJ, 720, 976 


A. Koch et al. 

Danforth, C.W., Tilton, E.M., Shull, J.M., et al. 2014, ApJ, 
submitted (arXiv:1402.2655) 

Fardal, M. A., Babul, A., Guhathakurta, P., Gilbert, K. M., & 
Dodge, C. 2008, ApJ, 682, L33 

Fardal, M. A., Weinberg, M. D., Babul, A., et al. 2013, MNRAS, 
434, 2779 

Ferguson, A. N. M., Irwin, J. J., Ibata, R. A., Lewis, G. F., 
Tanvir, N. R. 2002, AJ, 124,1452 
Gilbert, K. M., Guhathakurta, P., Kollipara, P., et al. 2009, ApJ, 
705, 1275 

Giroux, M. L. &; Shull, J. M. 1997, AJ, 113, 1505 
Grillmair, C. J. 2009, ApJ, 693, 1118 

Green, J. C., Froning, C. S., Osterman, S., et al. 2012, ApJ, 744, 

Haardt, F. & Madau, P. 2012, ApJ, 746, 125 
Helmi, A., & White, S. D. M. 1999, MNRAS, 307, 495 
Hendricks, B., Koch, A., Walker, M., et al. 2014, A&A, 572, AA82 
Huo, Z.-Y., Liu, X.-W., Xiang, M.-S., et al. 2013, AJ, 145, 159 
Ibata, R. A., Gilmore, G., Irwin, M. J. 1994, Nature, 370, 194 
Ibata, R., Irwin, M., Lewis, G., Ferguson, A. M. N., Tanvir, N. 
2001, Nature, 412, 49 

Ibata, R., Chapman, S., Ferguson, A. M. N., et al. 2004, MNRAS, 
351, 117 

Ibata, R., Martin, N. F., Irwin, M., et al. 2007, ApJ, 671, 1591 
Ibata, R. A., Lewis, G. F., McConnachie, A. W., et al. 2014, ApJ, 
780, 128 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, 
A&cA, 440, 775 

Keeney, B. A., Danforth, C. W., Stocke, J. T., Prance, K., & 
Green, J. C. 2012, PASP, 124, 830 
Keeney, B. A., Stocke, J. T., Rosenberg, J. L., et al. 2013, ApJ, 
765, 27 

Koch, A., Rich, R. M., Reitzel, D. B., et al. 2008, ApJ, 689, 958 
Layden, A. C., & Sarajedini, A. 2000, AJ, 119, 1760 
Lehner, N., Howk, J. C., Thom, C., et al. 2012, MNRAS, 424, 

Lehner, N., Howk, J. C., & Wakker, B. P. 2015, ApJ, 804, 79 
Lewis, G. F., Braun, R., McConnachie, A. W., et al. 2013, ApJ, 
763, 4 

Mackey, A. D., Huxor, A. P., Ferguson, A. M. N., et al. 2010, 

ApJ, 717, Lll 

Martin, C., Carlin, J. L., Newberg, H. J. & Grillmair, C. 2013, 
ApJ, 765, L39 

Mathewson, D. S., Cleary, M. N., Sc Murray, J. D. 1974, ApJ, 

190, 291 

Mayer, L., Mastropietro, C., Wadsley, J., Stadel, J., Sz Moore, B. 
2006, MNRAS, 369, 1021 

McConnachie, A. W., Irwin, M. J., Ibata, R. A., et al. 2003, 
MNRAS, 343,1335 

McConnachie, A. W., Irwin, M.J., Ibata, R.A., et al. 2009, 
Nature, 461, 66 

Mori, M., Sz Rich, R. M. 2008, ApJ, 674, L77 
Mould, J., &; Kristian, J. 1986, ApJ, 305, 591 
Newberg, H. J., Yanny, B., &; Willett, B. A., 2009, ApJ, 700, L61 
Osterman, S., Green, J., Froning, C., et al. 2011, Ap&SS, 335, 257 
Putman, M. E., de Heij, V., Staveley-Smith, L., et al. 2002, AJ, 
123, 873 

Rao, S. M., Sardane, G., Turnshek, D. A., et al. 2013, MNRAS, 
432, 866 

Reitzel, D. B., & Guhathakurta, P. 2002, AJ, 124, 234 
Rich, R. M., Mighell, K. J., Freedman, W. L., Sz Neill, J. D. 1996, 
AJ, 111, 768 

Richardson, J. C., Ferguson, A. M. N., Johnson, R. A., et al. 

2008, AJ, 135, 1998 

Richardson, J. C., Irwin, M.J., McConnachie, A.W., et al. 2011, 
ApJ, 732, 76 

Richter, P., Westmeier, T., Briins, C. 2005, A&A, 442, L49 
Richter, P., Charlton, J. C., Fangano, A. P. M., Bekhti, N. B., & 
Masiero, J. R. 2009, ApJ, 695, 1631 
Savage, B. D., &; Sembach, K. R. 1991, ApJ, 379, 245 
Searle, L. & Zinn, R. 1978, ApJ, 225, 357 
Sembach, K. R., &; Savage, B. D. 1992, ApJS, 83, 147 
Sembach, K. R., Wakker, B. P., Savage, B. D., et al. 2003, ApJS, 
146, 165 

Shull, J. M., Stevans, M., Danforth, C., et al. 2011, ApJ, 739, 105 
Stocke, J. T., Danforth, C. W., Shull, J. M., Penton, S. V., & 
Giroux, M. L. 2007, ApJ, 671, 146 
Thilker, D. A., Braun, R., Walterbos, R. A. M., et al. 2004, ApJ, 
601, L39 

Wakker, B. P., Sz van Woerden, H. 1991, A&A, 250, 509 
Wakker, B. P., Sz van Woerden, H. 1997, ARA&A, 35, 217 
Yanny, B., et al. 2003, ApJ, 588, 824 

York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 
120, 1579